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Abstract 

Hadron spectra and elliptic flow in relativistic heavy-ion collisions are studied in event-by-event 
(3+l)D ideal hydrodynamic simulations with fluctuating initial conditions given by the AMPT 
Monte Carlo model. Both the coherent soft gluon production from wounded nucleons and in- 
coherent mini-jet production in semi-hard parton scatterings are considered. These partons take 
part in parton cascade and are assumed to be locally thermalized. They provide the fluctuating 
initial conditions for hydrodynamic simulations with a Gaussian smearing. Effects of both trans- 
verse and longitudinal fluctuations are studied. The initial fluctuations along rapidity direction 
lead to expanding hot spots in longitudinal direction, which will reduce elliptic flow and the yield 
of particles at high transverse momentum. The intrinsic correlation introduced in the fluctuating 
initial conditions is also found to influence the di-hadron correlation of the final hadrons. 



1. Fluctuating initial conditions from AMPT 

Fluctuation in the initial condition for hydrodynamic study of heavy-ion collisions have been 
found necessary for the description of the anisotropic flow of final hadrons. We have carried 
out event-by-event (3+l)D ideal hydrodynamic simulations yj] with the local initial energy- 
momentum tensor in each hydrodynamic cell given by AMPT [2] model. AMPT uses HIJING 
Monte Carlo Model yj |4J for the initial parton production from hard scatterings and excited 
strings from soft interactions. The number of excited strings in each event is equal to participant 
nucleons. The number of mini-jets per binary nucleon-nucleon collisions follows a Poisson dis- 
tribution with the average value given by mini-jet cross section, which depends on the colliding 
energy and the impact parameter through an impact-parameter dependent parton shadowing |3J] 
in a nucleus. HIJING uses the Glauber model to determine the number of participants and binary 
nucleon-nucleon collisions with the Wood-Saxon nuclear distribution. Both soft partons from 
string melting and mini-jet partons will take part in parton cascade for a short period. In this 
model, the total local energy-momentum density of partons and its fluctuations will be deter- 
mined by the number of participants, binary nucleon-nucleon collisions, number of mini-jets per 
nucleon-nucleon collision, the fragmentation of excited strings, and the parton cascading. 

The 4-momenta and space coordinates of partons from the AMPT model according to the 
above description will be used to calculate the local energy-momentum tensor for the initial 



Preprint submitted to Nuclear Physics A 



November 8, 2012 



conditions of our event-by-event (3+l)D hydrodynamic simulations. Its value in each grid cell 
is approximated by a Gaussian distribution in invariant-time coordinates, 
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where pj = m iT cosh(T, - T] is ), pj = p ix , p] = p iy and p\ = m iT sinh(F/ - T} is )/T for parton i, 
which runs over all partons produced in the AMPT model simulations. We have chosen o> = 0.6 
fm, cr ;; = 0.6 in our calculations. The transverse mass mj, rapidity Y and spatial rapidity rj s are 
calculated from the parton's 4-momenta and spatial coordinates. Note that the Bjorken scaling 
assumption Y = rj s is not used here because of early parton cascade before the initial time and 
the uncertainty principle applied to the initial formation time in AMPT. We used K = 1 .45 and 
to = 0.4 fm in Au+Au collisions at V^aw = 200 GeV in the following calculations which gives 
good description of the multiplicity density of produced hadrons at central rapidity. 

2. The effects of longitudinal fluctuations on p T spectra and elliptic flow 

To investigate the effect of longitudinal fluctuations on hadron spectra and elliptic flow, we 
compare the (3+l)D hydrodynamic calculations using the full AMPT initial conditions with that 
using a tube-like smooth initial longitudinal distribution. In the event-by-event tube-like initial 
condition, we take the energy density and transverse flow velocity from AMPT results in the 
central rapidity region, and assume the fluctuation to be the same along longitudinal direction 
with an envelope distribution function, 



in rapidity. The length of the plateau rjo in the center rapidity region and the width of the Gaussian 
fall-off cr w at large rapidity are adjusted to fit the multiplicity distribution of charged hadrons. In 
another comparison, we also calculate the hadron spectra and elliptic flow for one shot tube 
initial condition, which is averaged over many AMPT initial conditions but rotated by an angle 
in each event to a common participant plane. We refer to this as one-shot-tube initial condition 
since in longitudinal direction the energy density and flow velocity also have a smooth tube-like 
distribution as defined in Eq. [2] 

Shown in Fig. [T]are the transverse momentum spectra (left panel) and elliptic flow (right 
panel) for n + in 30 - 40% semi-central Au+Au collisions at the RHIC energy yfsNN — 200 
GeV with the full AMPT initial conditions (solid lines), event-by-event tube initial condition 
(dot-dashed lines) and one-shot-tube initial conditions (dashed). The event-by -event fluctuations 
in the tube-like AMPT initial conditions significantly reduce the elliptic flow of final hadrons 
respect to the event plane as compared to the one-shot-tube initial conditions. The yield of high 
pj particles from event-by-event tube initial conditions on the other hand is increased due to the 
fluctuations (both in energy density and flow velocity) in the transverse direction as compared to 
the spectra from one-shot-tube initial conditions. Similar results were found in (2+l)D 15[] and 
(3+l)D [6] hydrodynamic simulations. This can be understood as the fast isotropic expansion of 
each hot spot in the transverse direction at the the early time of the hadrodynamic evolution. 

Fluctuations in the longitudinal direction in the full AMPT initial conditions also lead to 
hot spots in the longitudinal direction. Expansion of such longitudinal hot spots will dissipate 
more transverse energy in the longitudinal direction. This in turn decreases the elliptic flow at 
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Figure 1 : (Color online) The pj spectra (left panel) and elliptic flow (right panel) for charged pion given by hydrodynamic 
simulations at 30-40% semi-central Au+Au collisions for full AMPT initial conditions (solid lines), event-by-event tube 
initial conditions (dot-dashed lines) and one-shot-tube initial conditions (dashed lines). 



large pj as compared to the event-by-event tube initial condition. The yield of high pj particles 
are reduced due to the expansion of longitudinal hot spots from full AMPT initial conditions as 
compared to event-by-event tube-like AMPT initial conditions. 

Since the anisotropic flow is used to extract the transport coefficients (such as shear viscos- 
ity) from comparisons between experimental data and viscous hydrodynamics, the inclusion of 
longitudinal fluctuations in the hydrodynamic calculation will be necessary for more quantitative 
studies. 



3. The effects of longitudinal fluctuations on di-hadron correlation 

To study the effect of longitudinal fluctuations on the di-hadron correlation, we compare the 
di-hadron correlation for charged hadrons in (3+l)D hydrodynamic simulations with the full 
AMPT initial conditions with the event-by-event tube-like AMPT initial conditions. The two 
dimensional di-hadron correlation is defined as, 

CUiArj, A0) = <%Wfefe) (3) 

< N,(l]i, (pl)N a (l]2,(l)2) >mixed 

where N t and N a are the number of trigger and associate particles in two bins centered at {rj\, <p{) 
and (t]2, (pi)- Particles with the same rapidity A/7 = |?7i— 772I and azimuthal separation A<p = \<pi-cp2\ 
are averaged to increase statistics. The range of Arj and A(f> are extrapolated in the plot by 
symmetries of the correlation function. Correlation of particles from the same event is divided 
by that from mixed events to take into account of the finite rapidity window. In calculating 
the correlation for mixed events in the denominator, we rotate each event by a random angle in 
azimuthal direction to remove the background correlation coming from the same reaction plane 
used in AMPT model. 

We show in Fig. [2] the di-hadron correlation for charged hadrons in (3+l)D hydrodynamic 
simulation from full AMPT initial conditions and event-by-event tube-like initial conditions. The 
di-hadron correlation is flat in rapidity direction for tube-like initial conditions as compared to the 
full AMPT initial condition. This proves that the hydrodynamic simulation is not able to produce 
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(a)Full AMPT initial condition (b)Tube initial condition 




Figure 2: (Color online) Di-hadron correlations for charged hadrons calculated from (3+1) hydrodynamics for 30 - 40% 
semi-central Au+Au collisions at RHIC energy i/s^n = 200 GeV, from (a) full AMPT initial conditions and (b) event- 
by-event tube-like AMPT initial conditions. The transverse momentum cut-off for both trigger and associate particles 
are 2 < p T < 3 GeV/c. 

a near-side peak for low pj particles if no fluctuations and intrinsic correlations in longitudinal 
direction are introduced in the beginning of the hydrodynamic evolution. 

4. Summary 

Using fluctuating initial conditions from the AMPT Monte Carlo model for event-by-event 
(3 + 1 )D ideal hydrodynamic simulations, we introduce both fluctuations and intrinsic correlations 
in rapidity direction in the initial condition. Effects of longitudinal fluctuations (both in energy 
density and flow velocity) on hadron spectra, elliptic flow and di-hadron correlation are studied 
in hydrodynamic simulations from full AMPT initial condition, event-by-event tube-like initial 
condition and one-shot-tube initial condition. We showed that fluctuations in rapidity direction 
reduce noticeably the elliptic flow at large transverse momentum, which brings more constraint 
on extracting rj/s from comparison between the experimental data with viscous hydrodynamic 
calculations. The initial longitudinal fluctuations with intrinsic correlation built in is found to be 
important to understand the near-side peak of the di-hadron correlations for final particles at low 
p T . 
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